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ABSTRACT 

We present a new theoretical calculation of the contribution to the extragalactic gamma-ray background ra- 
diation (EGRB) from star-forming galaxies, based on a state-of-the-art model of hierarchical galaxy formation 
that is in quantitative agreement with a variety of observations of local and high-redshift galaxies. Gamma-ray 
luminosity (L 7 ) and spectrum of galaxies are related to star formation rate (i/f), gas mass (Af gas ), and star for- 
mation mode (quiescent or starburst) of model galaxies using latest observed data of nearby galaxies. We try 
the two limiting cases about gamma-ray production: the escape limit (L 7 oc ipM sa $) and the calorimetric limit 
(L 1 oc ijf), and our standard model predicts 7 and 4% contribution from star-forming galaxies to the total EGRB 
flux (including bright resolved sources) recently reported by the Fermi Gamma-Ray Space Telescope. System- 
atic uncertainties do not allow us to determine the EGRB flux better than by a factor of ~ 2. The predicted 
number of nearby galaxies detectable by Fermi is consistent with the observation. Intergalactic absorption by 
pair-production attenuates the EGRB flux only by a modest factor of ~1.3 at the highest Fermi energy band, 
and the reprocessed cascade emission does not significantly alter EGRB at lower photon energies. The sum 
of the known contributions from AGNs and star-forming galaxies can explain a large part of EGRB, with a 
remarkable agreement between the predicted model spectrum and observation. 
Subject headings: cosmic rays — diffuse radiation — galaxies: evolution — gamma rays : theory 



1. INTRODUCTION 

The existence of the extragalactic diffuse gamma-ray back- 
grou nd (EGRB) has been revealed first by the SAS-2 satel- 
lite dFichteletal.ll 19771: iFichtel et al]|1978l) . Better determi- 
nations of the flux and spectrum of EGRB were achieved by 
the EGRE T detector on board the Compton Gamma-R ay Ob- 
servatory (ISreekumar et all 1 19981: IStronget afl l2004a). The 
most reliable measurement of EGRB has very recently been 
reported b ased on the data of t he Fermi Gamma-Ray Space 
Telescope (lAtwood et al.l 120091) . and the EGRB spectrum is 
well described by a single power-law with a photon index of 
2.41 ± 0.05 and the photon flux of about (1.03 ± 0.17) X 
10~ 5 photons cm~ 2 s _1 sr _1 above 100 MeV dAbdo et al.l 

i2o!oa . 

The origin of EGRB has been discussed for a long time 
and various sources have been discussed as possible contrib- 
utors to EGRB, such as active galactic nuclei (AGNs, espe- 
cially blazars), galaxy clusters, intergalactic shocks produced 
by structure fo rmation, and dark matter annihilation (see, e.g., 
Dermer 2007 for a review). 

Almost all of the known extragalactic gamma-ray sources 
are blazars, and their contribution to the EGRB has been in- 
tensively studied (e.g., Inoue & Totani 2009; Venters 2010; 
Abdo et al. 2010d, and references therein). However, 
star-forming galaxies should also be gamma-ray emitters 
by cosmic-ray interactions with interste llar medium (ISM) 
and interste llar radiation field (ISRF) ([Strong et al.l l2000t 
Strong et al. 2004b), and there must be non-zero contribu- 
tion to EGRB. This is obvious because we know that the 
Galactic disk is a strong source of diffuse gamma-rays, 
and gamma-rays from Large Magellanic Cloud (LMC) have 
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been detected by EGRET (ISreekumar et al .1119921) . Further- 
more, gamma-ray emission in Ge V-TeV from Small Mag- 
ellanic Cloud (SMC, lAbdo et all 1201 Pel) and two nearby 
starburst galaxies, M82 and NGC 253, have recently 
been discovered by H.E.S.S (lAcero et all 120091). VERI- 
TAS ([VERITAS Collaborationll2009l) . and Fermi dAbdo et al.l 
l2010al) . The purpose of this paper is to present a new estimate 
of the contribution from star- forming galaxies to EGRB. 

There are a number of previous studies on this issue (Strong 
et al. 1976; Lichti et al. 1978; Dar & Shaviv 1995; Pavli- 
dou & Fields 2002; Thompson et al. 2007; Bhattacharya & 
Sreekumar 2009; Ando & Pavlidou 2009; Lacki et al. 2010; 
Fields et al. 2010), and the estimates of the contribution to 
EGRB ranges -10-50 ( %fl Most of these studies calculated 
gamma-ray luminosity (L 7 ) from star formation rate (SFR) 
of galaxies, because the cosmic -ray energy input is expected 
to be proportional to SFR. However, if escape of cosmic-rays 
from galaxies is signi fican t, SFR cannot be used as a reliable 



indicator of L 7 (see 



in more details). Several studies 



used infrared luminosity of galaxies as SFR indicators, but IR 
luminosity (especially at far-IR in early-type galaxies) is not 
a perfect SFR estimator because a significant amount of dust 
can be heated by ISRF from relatively old stars (e.g., Salim et 
al. 2009 and references therein). 

Amount of interstellar gas should also be important to de- 
termine L 7 , because the degree of cosmic -ray escape is af- 

1 Here, the contribution is against the total (physical) extragalactic back- 
ground photon flux (> 100 MeV) including bright resolved sources. It is often 
calculated against the unresolved component of EGRB, but it depends on the 
flux sensitivity of a particular detector. Throughout the paper, "the contribu- 
tion to EGRB" is defined against the total EGRB flux, which we estimated 
from the resolved and unresolved components of the Fermi data (Abdo et 
al. 2010d). We used a photon index of 2.4 for resolved blazars (Abdo et al. 
2010d) to extrapolate the resolved component from 200 to 100 MeV. 
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fected by the target amount. In fact, recent observations indi- 
cate that gamma-ray luminosity is nicely correlated with the 
product of SFR and gas mass in galaxies(Abdo et al. 2010a). 
Pavlidou & Fields (2002) included gas mass in the theoretical 
prediction of EGRB using the data of the cosmic star forma- 
tion history (CSFH). Galaxies in the universe is considered 
as a closed box containing stellar and gas mass in present- 
day galaxies, and CSFH is used to solve the time evolution 
of relative fractions of stars and gas. However this approach 
likely overestimates the gas mass contributing to gamma-ray 
production at high redshifts, because most of the present-day 
stellar mass is in the form of gas and assumed to contribute 
to gamma-ray emission. In reality, only the gas in collapsed 
dark halos can contribute to cosmic -ray interactions in galax- 
ies, which is a small fraction at high redshifts according to the 
structure formation theory. Recently, Fields et al. (2010) in- 
corporated gas mass in galaxies using the Schmidt-Kennicutt 
relation in their calculation of EGRB. In this case one must 
know galaxy size to relate SFR and Af gas , and a single value 
of galaxy size at each redshift was assumed as a simple model. 

Here we present a new calculation of the contribution from 
star-forming galaxies in EGRB using a state-of-the-art the- 
oretical model of galaxy formation in the framework of hi- 
erarchical galaxy formation, which is in quantitative agree- 
ment with a variety of observations at high redshifts as well 
as the local universe. An advantage of our approach com- 
pared with previous studies is that we can calculate both SFR 
and gas mass of individual galaxies at various redshits. Fur- 
thermore, we utilize the information of the gamma-ray spec- 
tra recently observed for nearby starbursts, in addition to the 
standard spectrum of the Galactic diffuse emission, to predict 
the EGRB spectrum based on the galaxy formation model in- 
cluding both quiescent and starburst galaxy populations. 

The paper is organized as follows. In the next section we 
describe our model calculations. We present the results on 
EGRB and statistics about the number of nearby galaxies de- 
tectable by Fermi in £|3] Af ter discussion on the uncertain- 
ties in our calcu lation ( £|4.U and implications for the origin 
of EGRB ( JJ4.21 . conclusions will be presented in $5] I n this 
work, cosmological parameters of Hq = 0.3, fl\ = 0.7, and 
Hq = 70 Mpc -1 km s _1 are adopted. 

2. FORMULATIONS 

2.1. The Theoretical Model of Hierarchical Galaxy 
Formation 

We use a mock numerical galaxy catalog produced 
by one of the latest so-called semi-analytic models 

i SAMs) of hierarchic al galaxy formation [the Mit aka model 
Nagashima & Yoshii 2004; N agashima et al.ll2005l) l. In gen- 
eral, SAMs compute merging history of dark matter (DM) ha- 
los based on the standard structure formation theory driven 
by cold dark matter, and include several important phys- 
ical processes related to the evolution of baryons in DM 
halos such as radiative gas cooling, star formation, super- 
nova feedback, galaxy merger, stellar population synthesis, 
chemical evolution, and extinction by inte r stellar dust (e.g., 
Kauffmann et al. 1993; Cole et al. 1994; Nagashima et al. 
1999; Somerville & Primack 1999; Baugh et al. 2005). In 
each time step of the halo merging history, two discrete modes 
of star formation are considered in the model, i.e., starburst 
and quiescent. When galaxies experience a major galaxy 
merger, intensive burst of star formation is assumed to occur 
(the starburst mode), and all of the available cold gas is con- 



verted into stars and hot gas within a short time. Otherwise 
star formation proceeds at a modest rate determined by cold 
gas amount and dynamical time scales (the quiescent mode). 

The Mitaka model can quantitatively reproduce a wide va- 
riety of observed characteristics of local galaxies, including 
luminosity functions (LFs) and scaling relations among vari- 
ous observables such as magnitude, colors, surface brightness, 
size, gas mass-to-light ratio, and metallicity (Nagashima & 
Yoshii 2004). Moreover, it can also reproduce well the cosmic 
star formation history (Nagashima & Yoshii 2004), the rest- 
frame ultraviolet (UV ) continu um LF of Lyman break galax- 
ies at z — 4 and 5 (Kashikawaet al. 2006), and all of the 
available observations for the high-z Lya emitters (Lya and 
UV continuum LFs, and Lya equivalent width distributions) 
dKobavashi et al.ll2007L[2010h . 

In the version of the Mitaka model used here, about 100 
Monte-Carlo realizations of merger histories are produced for 
DM halos in a DM mass interval of 0. 1 dex at each output 
redshift ranging z = 0-10. The model calculates DM ha- 
los with velocity dispersions larger than 30 km/s, and this al- 
lows us to make a reliable prediction in the stellar mass range 
of > 10 6 Af Q , roughly corresponding to the absolute B band 
magnitude of Mb < — 10 mag. 

2.2. Modeling Gamma-Ray Emission from Galaxies 

To calculate the EGRB flux from normal galaxies, we need 
the gamma-ray LF <f)(L 1 ,z) (comoving number density of 
galaxies per unit gamma-ray luminosity) as a function of red- 
shift. Here, L 1 is defined as gamma-ray luminosity in the rest- 
frame photon energy range of 0.1-5 GeV. The LF </> is calcu- 
lated from the mock galaxy catalog by modeling gamma-ray 
luminosity of each galaxy as follows. 

Diffuse gamma-ray radiation in galaxies is produced by the 
interactions between the cosmic -rays and ISM. If the cosmic- 
ray losses are dominated by escape, the gamma-ray luminos- 
ity is set by equilibrium between cosmic ray production and 
escape. It is then expected that gamma-ray luminosity de- 
pends on SFR (if), an indicator of cosmic ray production rate) 
and mass of interstellar gas (M gas , the amount of target for 
cosmic rays to interact before escape), i.e., L 1 oc ipM gas . 
We call it "the escape limit". However, if the cosmic -ray en- 
ergy losses are dominated by inelastic collision with interstel- 
lar gas and almost all of the cosmic-ray energy is converted 
into gamma-rays, the gamma-ray luminosity is no longer de- 
pendent on the amount of gas, i.e., L 1 oc ip. We call it "the 
calorimetric limit". [See also Pavlidou & Fields (2001), Tor- 
res et al. (2004), and Thompson et al.(2007) for more discus- 
sion on this issue.] Since it is difficult to model the detailed 
escape and energy loss processes of cosmic-rays in galaxies, 
we consider both of these two simple limiting cases. 

To establish the link between L 1 and ip and/or Af gas , 
we collected the observed values of these quantities for the 
nearby galaxies of SMC, LMC, the Milky Way (MW), M82, 
and NGC 253. They are summarized in Table Q] For 
SMC and LMC , SFR s are estimated from Ha lu minosities 
(Kennicutt et alj 120081) . using the relationship of iKennicuttl 
(1998) assuming the Salpeter initial mass function (IMF) in 
0.1-100 Mq. Extinction by dust is not taken into acc ount , 
and the uncertainty about this will be discussed later ( 34.U . 
For M82 and NGC 253, we used SFR estimated from the 
I R luminosity from interstellar dust, using the relationship 
of lKennicuttl ([1998) assuming the same IMF as ab ove , since 
these galaxies has large extinction (Dale et al. 20071 and lPencd 
1980) and Ha is no longer a good indicator of SFR. To- 
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tal IR luminosi ty is calculated f rom the IRAS fluxes in the 
three bands (Sanders et al. 2003) using the formula given in 
iDale & Hel ou (2002). Since it is difficult to estimate the total 
IR luminosity of MW, we use the supernova rate estimated by 
a gamma-ray observation (Diehl et al. 2006) and converted it 
into SFR assuming the same IMF. The standard value of 8 Mq 
is adopted for the mass threshold of core-collapse supernova 
explosions. 

The gas mass is estimated by the total of H2 and 
H I, where the former and the latter are measured by 
CO and 21 -cm line observations, respectively. The refer- 
ences are shown in Table [T] There are various estima- 
tions for H I mass of SMC an d LMC by differe nt au- 
thors (Huchtmeier&Richter 1989, Stanimirovi c et al.lfl999l 



Brims et al. 2005 for SMC, and Huchtmeier & Richter 1989, 
Luks & Rohlfs 1992, Westerlund 1997, Stavelev-Smith et al. 
2003, Brims et al. 2005 for LMC). Here we adopted the lat- 



est value by Bruns et al. (2005), and the standard deviations 
of the H I mass measurements among difference papers were 
added to the errors as a quadratic sum, to take into account the 
uncertainty (TableQ}. 

The collected data are plotted in Figure Q] It is clearly seen 
that there is a good correlation both for the L^-ipM gas (top 
panel) and L^-ip relations (bottom panel). We fitted these re- 
lations by a single power-law function to the data points of 
SMC, LMC, M82, and NGC 253, and the results are 



L 7 = (0.28 ±0.07) 



</, 



M ff 



M Q yr- 



10 9 Mr 



0.86±0.06 



[lO^erg/s] (1) 



and 



L 7 = (0.46 ±0.12) 



i> 



M Q yr- 



1.2±0.09 



[10 39 erg/s] 



(2) 



where the errors are statistical la. In the fitting calculation 
we adopted the effective variance method which takes into 
account the errors of both coordinates (Orear 1981). Since 
SFR of MW has a large uncertainty, we did not use the data 
of MW in the fitting, although it is consistent with the corre- 
lation. The exact proportionality between L 7 and ipM gas or 
L 1 and ip are not necessarily expected, because it depends on 
propagation and degree of confinement of cosmic rays within 
a galaxy. Therefore we use these power-law relation to cal- 
culate gamma-ray luminosity. The cases of the exact propor- 
tionality, i.e., L 7 oc (ipM gas ) 1 and L 1 oc (ip) 1 will also be 
discussed later (£ ]4.1I ). 



It should be noted that M gas in the Mitaka model has been 
compared with observations of local galaxies as a function of 
B band galaxy luminosity (Nagashima & Yoshii 2004; Na- 
gashima et al. 2005), showing a good agreement. This gives a 
support to apply A/ gas of the Mitaka model to the L 7 -^M gas 
relation determined by observations, though the model has not 
yet been tested against M gas at high redshifts because only 
very few observations are available. As for SFR, the adopted 
IMF in the Mitaka model (the Salpeter in 0.1-60 M Q ) is 
slightly different from that (in 0.1-100 Mq) assumed in con- 
verting the observed IR luminosity to SFR in Fig. [T] Though 
the value of SFR is hardly affected by the upper bound of 
stellar mass in IMF, it may affect the conversion from IR lu- 
minosity to SFR. However we confirmed, using the PEGASE 
stellar library (Fioc & Rocca-Volmerange 1997), that the dif- 



ference of the IMF mass ranges hardly affects the total lumi- 
nosity from young stars because stars heavier than 60 Mq are 
scarce and have short lifetimes. It is negligible at the stel- 
lar population ages larger than 10 65 yrs, and at most 20% at 
ages younger than that. Therefore we can apply SFRs in the 
Mitaka model to the fitting formula eqs. (Q} and Q, 

2.3. Gamma-Ray Spectra of Galaxies 

We also need to model the gamma-ray spectrum of galax- 
ies to predict the EGRB spectrum. The observed gamma-ray 
spectra of MW, M82, and NGC 253 are shown in Figure |2] 
We simply assume that quiescent galaxies have the same spec- 
tral shape as that of MW. For the MW spectrum we used the 
GALPROP fit to the Fermi spectral data of the Galactic dif- 
fuse gamma-rays dAbdo et al.l2010d) . The observed spectrum 
of SMC and LMC are similar to that of MW. 

For starburst galaxies, we expect some change of gamma- 
ray spectrum by the following physical reason. Because of 
the high density of star forming regions, it is plausible that 
starburst galaxies are "calorimetric" and have harder spectrum 
than quiescent galaxies because of less significant espace of 
high energy cosmic rays. Therefore we apply two modelings 
for the starburst spectra: (1) starbursts have the same spec- 
trum as MW, and (2) the MW-like spectrum at low energy 
range and a hard power-law tail at high energy range. We 
fitted the spectrum of M82 and NGC 253 by the MW spec- 
trum at energy less than 10 GeV, and with power law above 
the 10 GeV. Normalization of the power law component was 
determined so that it smoothly connects to the MW spectrum 
at 10 GeV. (Therefore the free parameters of the fit are two: 
the normalization of MW component and the photon index 
of high-energy power-law component.) As a result of the fit- 
ting, we find that M82 and NGC 253 can be fit with almost 
the same photon index of 2.2, and hence we adopt this spec- 
trum when we try a different spectrum from MW for starburst 
galaxies^. These model spectra are shown in Figure [2] 

2.4. The EGRB Calculation 

Given the modeling of gamma-ray luminosity of galax- 
ies described above, the EGRB flux and spectrum can be 
calculated by integrating the gamma-ray luminosity function 
4>{L 11 z) over redshift. To see the redshift evolution of the 
gamma-ray emissivity, we show the redshift evolution of p^ 
(comoving SFR density, proportional to gamma-ray luminos- 
ity density £ 7 in the calorimetric limit) as well as p^M (the 
comoving density of ipM gas , proportional to £ 7 in the escape 
limit) in Figure [3] This figure indicates that p^M decreases 
more rapidly toward higher redshift than p,/;, and this is likely 
because ip is roughly proportional to M gas and hence both 
the cosmic densities of ip and Af gas decrease to high redshift 
beyond z ~ 1. 

Given the gamma-ray luminosity function of galaxies 
<p{L~ n z), the EGRB spectrum (photon flux per unit photon 
energy per steradian) is expressed as 



d 2 F(e 1 ) _ c 
de 1 dfl 4n Ju 



dt 

dz 



dz I dLj4>(L 7l z) 




x {l + z) dL,[L^(l + z)e,} e _ Tj{z ^ ^ (3) 
where e 7 is photon energy for observers at z = 0, t the cosmic 

2 In the spectral fitting, we ignored the data points that have only upper 
limits. 
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Table 1 

Summary of L 7 , SFR, M gaa , and Distance of Nearby Galaxies 



10- 



Object 



L 7 
[10 39 erg/s] 



— smc — 0.011 ±0.001 
LMC 0.041 ± 0.007 
M82 13.0 ±5.0 

NGC 253 7.2 ± 4.7 

MW 3.2 ± 1.6 



SFR" 

[M Q /yr] 



0.037 ±0.011 

0.244 ±0.073 

16.3 ±2.4 

7.9 ±4.9 

2.6 ±1.5 



M g 
[1O H M ] 



0.45 ±0.11 
0.56 ±0.14 
4.9 ±0.58 
4.3 ±0.55 
7.0 ±1.0 



Distance 
[Mpc] 



0.061 ± 0.003 

0.049 ± 0.001 

3.6 ±0.3 

3.9 ±0.4 



a Gamma-ray luminosities are in 0.1-5 GeV, from Abdo et al. (2010e) for 
SMC, Abdo et al. 1 2010b) for LMC, and Abdo et al. 1 20 10a) for M82, NGC 
253 and MW. 

b See text ( £|2.2t for methods and references of SFR estimates. 

c Total gas mass (H I + Hg) are estimated from 21-cm and CO line 
fluxes. References are WeiB et al. 12005), Huchtmeier & Richter 1 1989) (for 
M82), Knudsenetal. (2007), Springob et al. (2005) (for NGC 253), and 
Boissier & Prantzos ( 1999) (for MW). For SMC and LMC, see text ( EH 
for methods and references of gas mass estimates. 

"■References are Hilditch et al. (2005) for SMC, Pietrzyhski et al. (2009) for 
LMC, Freedman et al. 1 1994) for M82, and Karachentsev et al. (2003) for 
NGC 253. 

time, and dt/dz can be calculated if the standard cosmolog- 
ical parameters are given. The rest-frame gamma-ray spec- 
trum of galaxies cLL 7 (L 7 , t 1:T )/dt 1:r can be calculated for a 
given L 7 by the spectral shapes as modeled in the previous 
section, where e 7 , r is a photon energy at the rest- frame of the 
source redshifts. We set the maximum redshift z max = 10 
corresponding to that of the galaxy formation model, but this 
parameter hardly affects the EGRB flux because the major- 
ity of the EGRB flux comes from galaxies at z < 2 (see S]3] 
below). 

High energy photons (> 20 GeV) are absorbed during the 
propagation in intergalactic space, by interaction with the cos- 
mic infrared ba ckground (CIB) photons producing electron- 
positron pairs (iSalamon & Sti cker 1998; Totani & Takeuchil 
|200H iKneiske et al.l I2004t iStecker et a l. 2006). The opti- 
cal depth r 7 (z, e 7 ) in equation (2) represents this interaction. 
Since the infrared emission from dust has not yet been in- 
cluded in the Mitaka model, r 7 cannot be calculated self- 
consistently based on the Mitak a mode l. Therefore we use 
the model of Totani & Takeuchi (2002) for r 7 in this work, 
which is based on a different galaxy evolution model. To ex- 
amine the effect of the uncertainty about r 7 , we also calcu- 
lated with differen t mode ls of r 7 by IKneiske et al. (2004) and 
Franceschi ni et al.l (120081) . and found that the change in the 
EGRB flux is not significant (negligible at < 30 GeV, and at 
most by a factor of 1.5 at the highest Fermi energy band). 

The created pairs would scatter up the cosmic microwave 
background photons to gamma-rays via inverse Compton 
mechanism (cascade emission) and they make some con- 
tributions to the EGRB flux at low er photon energies 
dAharonian et al.lll994t iFan et al.l 120041) . We calculate the 
contributio n from the cascade em ission using the same for- 
malism of Kneiske & Mannheim (2008). 

3. RESULTS 

3.1. EGRB Predictions Compared with Observation 

Figure @] shows the quiescent and starburst components of 
EGRB calculated by the formulations presented in the previ- 
ous section, for the two different modelings of the escape and 
calorimetric limits in different two panels . The data points ar e 
the Fermi measurements of EGRB flux dAbdo et al.ll2010cT) . 
The reported Fermi EGRB spectrum is the residual after the 



LMC 



SMC 




V[M Q yr" 



Figure 1. (Top) Correlation between the gamma-ray luminosity L 7 (0.1-5 
GeV) and the product of SFR and gas mass, ipM gas . The solid line represents 
the best fit power-law L 1 <x (i/>Af gas ) a with a taken as a free parameter. 
The dashed line represents a linear function fit with a fixed value of a = 1. 
(Bottom) The same as top panel, but for the correlation between L 7 and SFR. 



sources resolved by Fermi have been subtracted, and it is dif- 
ficult to compare a theoretical prediction to such data tak- 
ing into account complicated detection process and efficien- 
cies. Therefore we also plot the total Fermi EGRB flux in- 
cluding the cont ribution from resolved sources as reported by 
lAbdo et ail d2010d) . and compare our prediction including all 
sources in the universe with this total observed flux. For the 
starburst galaxies we plotted the two models using the two dif- 
fer ent s pectral shapes (MW or MW+power-law) as discussed 
in J12.3I All the predictions take into account the intergalactic 
absorption and the reprocessed cascade emission as discussed 
in the previous section. 

In Figure[5] the predictions of total EGRB flux from all star- 
forming galaxies (quiescent plus starburst) are plotted. To 
show the effect of intergalactic absorption, the EGRB spec- 
trum without taking into account the absorption is also plot- 
ted. The cascade component produced by the absorption is 
also shown. It can be seen that the intergalactic absorption at- 
tenuates the EGRB flux by a factor of about 1 .3 at the highest 
photon energy of Fermi. The fraction of the absorbed energy 
flux of EGRB is not a large amount compared with the total 
EGRB energy flux, and hence the cascade component does 
not make a significant contribuiton at the low energy bands 
(Note that Fig. [5]is made as a vF v plot). 

The EGRB flux from starburst galaxies is significantly en- 
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10 io z 

Energy [GeV] 

Figure 2. The gamma-ray spectrum E 2 dN/dE of MW, M82 and NGC 253, 
where E is gamma-ray energy and dN/dE the differential photon flux. The 
MW spectrum is diffuse background flux in units shown on the right ordi- 
nate, and the others are flux from a source in units shown on the left ordinate. 
For the purpose of presentation, the spectrum and the data points of M82 
are multiplied by 100. The MW data (open squares) are the Galactic diffuse 
gamma-rays averaged over \b\ > 10° I Abdo et al. 2010c). The Fermi data of 
M82 (fi lled circles) and NGC 253 (open circles) are obtained by Abdo et al. 
1201 Oat). The data in the very high-energy range of M82 (filled triangles, 
VERITAS Collaboration 2009) and NGC 253 (open triangles, lAcero et all 
2009) are also plotted. The solid curves are the model spectra used to calcu- 
late EGRB (see text for details). 
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Figure 3. The evolution of the comoving densities of SFR (p^, , solid 
and the product of SFR and gas mass (p^m> dashed line) in galaxies, 
culated by the model of hierarchical galaxy formation. See the left and 
ordinates for the units of p^ and p^j^r, respectively. 
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hanced at high energy range of e 7 > 10 GeV when the 
MW+power-law spectra is applied. However, the differ- 
ence is not significant in the total EGRB flux from all star- 
forming galaxies, provided that quiescent galaxies have simi- 
lar gamma-ray spectra to that of MW. It should also be noted 
that the classification of quiescent and starburst galaxies in the 
Mitaka model is determined to reproduce the optical/infrared 
observed properties of galaxies. Therefore the classification 
can be systematically different if we do it by the difference of 
gamma-ray spectra. This is an interesting question but diffi- 
cult to answer at the present time because of too few observed 
data of gamma-ray spectra of galaxies. 

The EGRB photon flux from all star-forming galaxies is 
9.9 xlO (escape limit) and 5.6 xlO (calorimetric limit) 



photons cm" 2 s" 1 sr" 1 (> 100 MeV), which are 7.0±1.8% 
and 3.9±1.0% of the EGRB flux observed by Fermi (includ- 
ing resolved sources), respectively. If we compare these re- 
sults with the recent Fermi measurement of the unresolved 
EGRB flux, the contribution from star-forming galaxies be- 
come 9.6±2.4% (escape limit) and 5.4±1.4% (calorimetric 
limit). Here the error is that coming from the statistical lcr 
error of the fitting formula of L 1 [eqs. ((TJ and (0]. 

As a result of calculation the EGRB flux from star-forming 
galaxies in the calorimetric limit becomes lower than the es- 
cape limit. This seems counter-intuitive, since more gamma- 
rays should be produced in the calorimetric limit if the 
cosmic -ray energy input is fixed. This is because we use the 
datapoints of SMC and LMC to derive eq. (0. These galax- 
ies are likely in the escape limit, making the slope of ip-L^ 
relation steeper. It would be an overestimate of EGRB if we 
calculate the "truly" calorimetic case, i.e., L 1 ex ip 1 with the 
proportionality constant for starburst galaxies. Rather, what 
we calculated here is the best empirical estimate of EGRB in 
a model in which L 1 is a function of tp, 

3.2. Redshift Distribution of EGRB Photons 

Figure [6] shows the cumulative redshift distribution of the 
EGRB photons from star-forming galaxies. It is clearly seen 
in the figure that more than half of the total EGRB flux comes 
from galaxies at z < 1, and more than 90% from z < 2 
both for the escape- and calorimetric-limit models. This is a 
reasonable resu lt given the evolution of p^, and p^m as dis- 
cussed in i )2.4l The distribution of photons from starburst 
galaxies extend to slightly larger redshifts than that from qui- 
escent galaxies, because of the stronger evolution of p^M at 
z ~ 0-1. Our result of the escape-limit can be compared with 
that of Ando & Pavlidou (2009, hereafter AP09). The two 
models are in rough agreement, though the median redshift 
of our model is slightly smaller than that of AP09. This may 
be because AP09 assumes that all the mass in the present-day 
stars contributes to gamma-ray production as gas in the early 
universe, while only a small fraction is in the form of cooled 
gas in our model, as discussed in Sjl] 

3.3. Detectability of Nearby Star-Forming Galaxies 

We also calculated the number of galaxies that are de- 
tectable by Fermi. The 1-yr sensitivity of Fe rmi is ~ 3 x 
10~ 9 photons cm _2 s _1 at e 7 > 100 MeV dAtwood et al.1 
2009). Our model predicts 4.6 and 2.3 galaxies above this 
flux in the escape- and calorimetric-limit models, respectively. 
These numbers are consistent with the actual number of four 
(SMC, LMC, M82, and NGC 253) or two (when SMC and 
LMC are excluded because it is within our Galaxy halo; see 
below). This result gives another support to the reliability of 
our EGRB prediction. However, it should be noted that in 
this estimate all galaxies in the Mitaka model are assumed 
to distribute randomly in space with a uniform mean den- 
sity, because the clustering information of galaxy distribution 
is not included in the semi-analytic models like the Mitaka 
model. The density fluctuation and clustering in the local 
group around MW should affect the above estimate, espe- 
cially for satellite galaxies like LMC. Not only the cluster- 
ing effect but also the realistic detection efficiency of Fermi 
sources around the detection limit must be properly taken into 
account in a more quantitative analysis, which is beyond the 
scope of this paper. In the future, Fermi sensitivity will reach 
- 1 x 10~ 9 photons cm" 2 s" 1 (> 100 MeV) by the 10- 
year observation, and we expect that Fermi will detect ~ 24 
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(escape limit) or 12 (calorimetric limit) nearby star-forming 
galaxies, though it is again subject to the uncertainty about 
clustring, and is probably optimistic assuming 100% detec- 
tion efficiency above that flux. 

Since the Fermi sensitivity becomes worse for extended 
sources than that for point sources, we checked the size of 
galaxies computed by the Mitaka model. There is a large 
scatter in the correlation between size and gamma-ray lumi- 
nosity, but the average angular size of model galaxies brighter 
than Fermi 1-yr sensitivity is ^0.05 deg, which is much 
small er than the Fermi a ngular resolution (~ 0.6 deg at 1 
GeV, lAtwood etaf] 12009). Therefore the effect of extended 
sources should not significantly affect the above estimate for 
the number of detectable sources. It should also be noted that 
the size of an active star-forming region in a galaxy can be 
much smaller than that of the stellar disk of the same galaxy, 
as inferred from observations of NGC 253 (Acero et al. 2009). 
Therefore gamma-ray emitting regions can be much smaller 
than the size of the model galaxies in the Mitaka model. 

4. DISCUSSION 

4. 1 . Uncertainties 

Here we examine the uncertainties about our estimate of 
the star -forming galaxy contribution to EGRB. As described 
in £12.21 we have calculated gamma-ray luminosity of galaxies 
by the relations L 7 oc (V^gas) ' 86 and L 1 oc (ip) 1 which 
are obtained by the fit to the data of SMC, LMC, M82, and 
NGC 253. We also tested the cases of L y oc (V'Afgas) 1 and 
L 1 oc (ip) 1 , and in this case, the best-fit relations to the data 
becomes 



L 7 = (0.36 ± 

.( V> 



M e 



^0 yr' 



-1 



To^r) l 1 " 39 ^/* < 4 > 



and 



i 7 = (0.28 ± 0.08) x ( — t — \ [10 39 erg/s]. (5) 
\Mq yr J 

The eq. Q is consistent with the estimate of the same re- 
lation by iPavlidou & Fieldsl (|2001) within the uncertainty of 
the conversion of MW supernova rate into SFR. We found 
that in these cases the EGRB flux becomes slightly differ- 
ent from our standard models: 9.3±2.0% (escape limit) and 
1.8± 0.5% (calorimetric limit) contribution to the observed 
total Fermi EGRB flux (or 12.8±2.8% and 2.5±0.7% against 
the unresolved component of EGRB reported by Fermi). The 
value ofl2.8±2.8% can be compared with another recent re- 
sult by Fields et al. (2010), based on more empirical calcu- 
lations about galaxy eovlution, and they are consistent with 
each other within the systematic uncertainties. 

Estimate of SFR is generally affected by the treatment of 
extinction by interstellar dust. The SFRs of LMC and SMC 
were estimated by their Ha luminosities without extinction 
correction. If a part of ionizing photons are absorbed by dust 
before hydrogen ionization, their energy is converted into IR 
luminosity. Therefore, the sum of SFRs estimated by Ha and 
IR luminosity gives a good upper bound on the total SFR. 
(It is an upper bound because non-ionizing photons also con- 
tribute to IR luminosity when they absorbed by dust.) The 
SFRs estimated by IR luminosities of SMC and LMC are 
0.0153±0.0046 and 0.146±0.04 M Q /yr, respectively, by the 



same meth od a s for SFR estimates of M 82 and NGC 253 de- 
scribed in £12.21 These are 40-60% of SFRs by Ha, indicating 
that these galaxies are not heavily obscured by dust. When we 
add the SFRs estimated by IR luminosity to the SFRs of SMC 
and LMC, the EGRB flux from star-forming galaxies changes 
slightly by a factor of 1.2. On the other hand, starburst galax- 
ies (M82, NGC 253) have large extinctions and hence SFRs 
estimated by IR l umin osity are much bigger than those by Ha, 
as mentioned in £12.21 

In the previous version of this paper put on the preprint 
server ( arXiv: 1005. 1390^1), we reported about 14% contribu- 
tion from star-forming galaxies to the total EGRB flux, which 
is about two times bigger than our new result. We exam- 
ined the reason for this, and found that it is mainly due to the 
smaller SFR of LMC used in the previous version, which was 
estimated only by IR luminosity. As mentioned above, SFR 
from Ha luminosity is about two times larger than that from 
IR for LMC ans SMC. (Adding SMC to the nearby sample is 
also a new point of the current version, but it does not signif- 
icantly affect the EGRB flux.) Therefore we believe that our 
new result is more appropriate, but it would be better to con- 
sider that the systematic uncertainty of the EGRB flux from 
star-forming galaxies is not smaller than a factor of two, given 
the small number of nearby galaxies with observed gamma- 
ray luminosities and uncertainties in estimates of SFR and gas 
mass. 

4.2. Implications for the Origin of EGRB 

Our result indicates that about 10% of EGRB can natu- 
rally be explained by star-forming galaxies. Here we discuss 
the origin of EGRB in general, also considering the larger 
contribution from AGNs. In Figs. [4] and [5] we plotted the 
EGRB model of AGNs by Inoue & Totani (2009, hereafter 
IT09). The IT09 model includes two populations of AGNs: 
one is blazars, and the other is non-blazars that are responsi- 
ble for the EGRB at 1-10 MeV. Here we plotted models of 
U03(g, 71) and Y = 3.5 (see IT09 for details). 

Many blazars have been detected by EGRET and Fermi, 
and it is obvious that they are the major contributor to EGRB. 
The blazar component of IT09 can account for 43% of the 
total EGRB flux, which is in good agreement with the latest 
observational estimate of Abdo et al. (2010d) based on the lat- 
est Fermi datqj. Therefore, adding the component from star- 
forming galaxies, more than 50% of EGRB can be explained 
by the sum of surely existing sources. The non-blazar AGN 
component of IT09 has not yet been confirmed by observa- 
tions, but it is introduced to explain the background radiation 
from hard X-ray to MeV region in a physically natural way 
(Inoue et al. 2008). If this component extends to ~ 100 MeV, 
the contribution to EGRB flux at > 100 MeV could be ~ 20%, 
and hence it is not difficult to explain > 70% of EGRB by the- 
oretically reasonable sources. It should also be noted that both 
the AGN component and star-forming galaxy component have 
remarkably similar spectra to the observed Fermi EGRB. 

The <30% residual of EGRB may not be explained simply 
by the sum of AGNs and star-forming galaxies, and there may 
be a room for different populations such as galaxy clusters or 
dark matter annihilation. However, there are still uncertainties 

3 Note that the number (16%) quoted by Abdo et al. (2010d) is the frac- 
tion of unresolved blazars against the unresolved diffuse EGRB observed by 
Fermi. If this number is converted to that of all blazars against the total (un- 
resolved + resolved components) EGRB flux, it increases to ~40% in good 
agreement with IT09. 
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Figure 4. The EGRB spectrum. The theoretical predictions for the contributions from quiescent and starburst galaxies are separately shown, for the escape-limit 
(left panel) and the calorimetric limit (right). For starburst galaxies, we plotted two cases corresponding to two different spectral shapes of MW and MW + 
power-law. For comparison, the EGRB models of blazars and all AGNs (blazar + non-blazar) by Inoue & Totani (2009) are also shown. Open circles show the 
total EGRB flux reported by Fermi ( Abdo et al. 2010c) including resolved sources, and the filled circles are the same but for the unresolved diffuse flux excluding 
resolved sources. The model curves should be compared with the open circle data. 
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Figure 5. The same as Fig. [4] but for the sums of each component are shown. The Models 1 and 2 are EGRB from all star-forming galaxies (quiescent plus 
starburst), with the two different spectra for starbursts. To show the effect of intergalactic absorption, the Model 1 without absorption (thin solid line) and the 
cascade component of the Model 1 are also shown. The total of star-forming galaxies plus AGNs is also shown using the Model 1 for star-forming galaxies. 



in the modeling of AGN gamma-ray luminosity function. The 
blazar luminosity function of IT09 was based on EGRET data, 
including uncertainties about detection and identification ef- 
ficiencies depending on the location in the sky. The recent 
observational estimate of the blazar contribution by Fermi 
data includes a large correction about detection efficiency near 
the sensitivity limit, which is generally dependent on models. 
Blazars are highly variable sources, and it should add another 
uncertainty in EGRB estimates. Therefore we conclude that a 
large part of EGRB can be explained by known or physically 
reasonable sources, and there is no strong evidence for exotic 
components like dark matter annihilation. 

5. CONCLUSIONS 

We have presented a new calculation of EGRB from 
cosmic-ray interactions in star-forming galaxies, based on 
a state-of-the-art galaxy formation model in the framework 



of hierarchical structure formation. The galaxy formation 
model is quantitatively consistent with various observations 
at high redshifts as well as the local universe. Gamma- 
ray luminosities of galaxies are calculated by star forma- 
tion rate and gas mass in model galaxies, with the relation 
L 7 oc (V>Mg as ) - 86 (the escape limit) or L 7 oc (^) 12 (the 
calorimetric limit) which are calibrated by the recent obser- 
vational data of nearby galaxies. The predicted number of 
nearby galaxies that are detectable by Fermi is consistent with 
the actual number observed so far. 

We found that star-forming galaxies make 4-7% contribu- 
tion to the total EGRB flux reported by Fermi in our stan- 
dard model. Combined with the contribution from blazars 
as estimated by the Fermi data, more than ^50% of EGRB 
can be accounted for. If the soft power-law tail of CXB is 
extending from MeV to GeV region as expected from the 
MeV background data and theoretical considerations of AGN 
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Figure 6. Cumulative redshift distribution of the EGRB photons (0.1—5 
GeV) from star-forming galaxies. Quiescent and starburst galaxies are sepa- 
rately shown by solid and dashed curves, respectively. 

accretion disks (Inoue et al. 2008), additional ^20 % can 
be explained. The combined spectrum by AGNs and star- 
forming galaxies is remarkably similar to the observed EGRB 
spectrum. It should be noted that there is no free parame- 
ters that can be tuned to fit the obeserved EGRB spectrum 
in the present model; the blazar EGRB spectrum of IT09 
is determined by the spectral templates of the blazar SED 
(spectral energy distribution) sequence, and that of the star- 
forming galaxy component in the present work by templates 
constructed by observed nearby galaxies. Therefore we con- 
clude that a large part (>70%) of EGRB can be explained 
by reasonable sources of AGNs and star-forming galaxies. 
Further examination is required to see whether the apparent 
<30% residual of EGRB is mainly a result of modeling un- 
certainties, experimental/observational uncertainties in deriv- 
ing the EGRB data, or significant contributions from com- 
pletely different source populations. Given the good spectral 
agreement of the Fermi data and AGNs/star- forming galaxies, 
the rest of EGRB is also likely dominated by astrophysical 
sources accelerating particles, even if a completely diffrent 
population is responsible for it. 

It is in contrast that a large part of the blazar component 
of EGRB will be resolved into discrete sources by the ul- 
timate Fermi sensitivity in the near future (IT09), while al- 
most all of the star-forming galaxy component will remain 
unresolved because of the faintness of individual sources. 
Therefore, any exotic contribution like dark matter annihi- 
lation cannot be probed directly under the level of the star- 
forming galaxies, i.e., ~ 10% of the total EGRB flux. An- 
other approach such as utilizing anisotropy would be required 
to sea rch for the signal under that level (see, e.g. Ando et al . 
2007: ISiegal-Gaskins & Pavlidoull2009l iHenslev et al.1 120091 
ICuocoetal.1120101) 

We assumed a simple empirical relation between gamma- 
ray luminosity, SFR, and gas mass of galaxies. Only two 
spectral templates were used in the calculation. A next step 
would be to construct a more physical model of gamma-ray 
luminosity and spectrum based on a larger number of phys- 
ical quantities (e.g., size and stellar radiation field in addi- 
tion to SFR and gas mass), taking into account propagation of 
cosmic -rays and production processes of gamma-rays. 
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